Positron annihilation in e"*"— H collision above the Positronium formation threshold 
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A long-standing problem to account for the electron-positron annihilation in positron Hydrogen scattering 
above the Positronium formation threshold has been resolved by the use of the three-body Faddeev formalism. 
The multichannel three-body theory for scattering states in presence of a complex absorbing potential is devel- 
oped in order to compute the direct e^e~ annihilation amplitude, the amplitude of Positronium formation and 
respective cross sections. A number of e^e~ direct annihilation cross sections and Positronium formation cross 
sections in the energy gap between Ps(ls) and H(n = 2) thresholds are reported for both the positron-Hydrogen 
incoming channel as well as the proton-Positronium incoming channel. 
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When colliding with the Hydrogen atom in the energy re- 
gion between Positronium Ps(ls) and Hydrogen H(n = 2) 
thresholds, the positron can annihilate in two possible chan- 
nels either in the elastic e+— H one or in the rearrangement 
p— Ps one. During reaiTangement the positron captures the 
electron from H and forms Positronium in the Ps(ls) state 
and then the positron-electron pair annihilates from this state. 
The lifetime of Ps(ls) is known to depend on the total spin 
(l|. This well known annihilation mechanism is of the two- 
body nature, i.e., once Positronium is formed the third par- 
ticle does not affect the annihilation process. The annihila- 
tion mechanism in the elastic e+ — H channel is more compli- 
cated and needs three-body footing. There was a number of 
attempts for the unified three-body treatment of the positron 
annihilation phenomenon and the Positronium formation, es- 
pecially in the vicinity of the Positronium formation threshold 
Q]^ 0> 01- Nevertheless as to the best of our knowledge, 
the self-consistent multichannel theory for the annihilation in 
positron-atom collisions when the reaiTangement channel is 
open does not exist to date. It is the aim with this paper to fill 
out this gap by combining the three-body multichannel scat- 
teringapproach basing on modified Faddeev equations (MFE) 
la], with the concept of a complex absorbing potential 
121 which allows to treat the annihilation in the framework of 
non-relativistic quantum mechanics. In the paper we restrict 
ourselves within the Ore gap energy region. The extension 
onto the case when more than one Hydrogen (and Positron- 
ium) state are participating is straightforward. 

Let us recapitulate those portions of the one channel ap- 
proach to the description of the annihilation phenomenon in 
positron Hydrogen scattering below the rearrangement thresh- 
old which are essential for a multichannel extension. The well 
known QED formula 121, M 

- ^rl{c/v)Z^^^ (1) 

expresses the spin averaged annihilation cross section ct" = 
1/4^ cr° + 3/4'' (7° in terms of the effective number of elec- 
trons ^eff contributing to the annihilation. The respective 
number ^'^'^^■^eff ^'^^ ^ ^P^'^ state. S is given by the integral 



''+'Z^^^= j dv,dv,\'S+^^\v,,v,)\^5{v,-v,). (2) 

Throughout the paper we use bold letters for vectors, e.g., 
r, not bold ones for their magnitudes, e.g., r = |r|, and we 
choose units such that h = 1. In formulae Q and (|2j the in- 
gredients has the following meaning, tq is the classical elec- 
tron radius, c is the velocity of light, v is the incident veloc- 
ity of positron, 2S+1yj,o jg jj^g p+ w scattering wave func- 
tion when only Coulomb interaction is taken into account, 
and ri and are position vectors of positron and electron 
in the coordinate frame associated with proton. The presence 
of the delta-function in (|2j manifests the fact that e+ and 
annihilate when they occupy the same position. The non- 
relativistic quantum mechanical description consistent with 
the above formulation can be based on the concept of a com- 
plex absorbing potential |2], |3]. The use of this potential in 
the Schrodinger equation leads to the particle loss. The cross 
section associated with this loss is actually the annihilation 
cross section. The absorbing potential is defined by 0] 

igW{Yi - ra) = -i(eVao)(''^+'c)a'5+3,5(ri - rg), (3) 

where e is the electron charge, ao is the Bohr radius, a is the 
fine structure constant, and spin dependent constants ^'^+^c 
are defined as = 27r and ''c = 8(7r^ — 9)/9. 

The Schrodinger equation for the positron Hydrogen wave 
function with outgoing asymptotic boundary conditions 

where is the three-body Coulomb Hamiltonian and pi is 
the momentum of incident positron, can be rewritten as 

= *pl + - H^r^w-^i, ■ (4) 

Here \I/p+ is the solution to the pure three-body Coulomb 
problem - = and = + iO. It is seen 

from that the asymptotics of , as ri approaches infin- 
ity, takes the form 

*+ (ri,r3) ^ ^A(r3){e'P^-""^ + [Fe - igFa]}. (5) 
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In Q, ^'(ra) is the Hydrogen ground state wave function. The 
elastic scattering ampUtude is entirely due to the 4*°^ con- 
tribution, and the annihilation amplitude Fa is due to the ab- 
sorbing potential. This amplitude is given by 

Fa{p\,Pi) = ^{^"^jW\^;,), (6) 

where m is the electron mass. In principle, the annihilation 
cross section cr" can be computed from the total cross section 
CT* as 

(7°- ^a* - j dfl \Fe\^, cr* = J dfl \Fe ~ igFa\^ , (7) 

but in the one channel case there is an alternative way which 
is based on the optical theorem 

at = — 3'm[Fe(pi,pi) - igFa{pi,pi)], 
Pi 

from which the annihilation cross section appears as 

47r 

cTa = g — 3m[-zFa(pi,pi)]. (8) 
Pi 

Now, the annihilation cross section Q and the formula (|3 for 
the effective number Z^ff follow immediately from (|6j and (|8} 
with replacing ^f^^ in (|6j by '1'°+ which, considering Q per- 
turbatively in view of jgl ^ 1, is the first order approximation 
forvI.+ ,/.e.,M'+ =*0t + O(.g). 

Let us now turn to the case when the Positronium forma- 
tion channel is open and thus the e+— H collision becomes 
truly multi arrangemental. In this case the integral defined by 
diverges and the cross section defined by formula Q is 
infinite. On the contrary, the exact amplitudes F^ and Fa de- 
fined by (|5}, <|6} are finite (as long as the exact solution to (|4} 
is assumed) due to the effect of the absorbing potential which 
makes the Positronium binding energy complex. Although the 
exact amplitudes are well defined, the matrix element (|6} can- 
not remedy the problem with the annihilation cross section. 
Indeed, the contribution of elastic and rearrangement chan- 
nels to the optical theorem for the three body system 1^ is 
additive 

47r 

— 3m[Fe(pi,pi) - igFa{pi,pi)] = an + a2i. (9) 
Pi 

Here an is the total (elastic plus annihilation) cross section 
in the e+— H channel and (72 1 is the total Positronium for- 
mation (and annihilation after Positronium formation) cross 
section. Since the contribution of the annihilation in the left 
hand side of (|9} is proportional to the matrix element (|6jl it is 
apparent that with this quantity the direct annihilation process 
and the process of the annihilation after Positronium forma- 
tion are inseparable. We believe that a variant of this effect of 
inseparability was observed in the numerical study of positron 
Hydrogen scattering in 

So, it is clear that the shortcoming of the outlined formalism 
is due to an inadequate treatment of multichannel scattering. 



The approach which is proven to define correctly the transi- 
tion amplitudes between all open channels in the three-body 
case is the Faddeev equation formalism |5] and its modifica- 
tion MFE for long-range interactions |6]. In this approach 
the wave function is split into the sum of components each 
of them corresponding to the given asymptotic arrangement. 
The solution of the set of equations for components (Fad- 
deev equations) satisfies the following requirements: i) the 
sum of components obeys the Schrodinger equation; ii) the 
asymptotics of the component associated with a given asymp- 
totic arrangement includes contributions from all open chan- 
nels in this arrangement and coincides with the restriction of 
the wave function on this arrangement |6]. This last property 
guaranties that the amplitudes computed from components are 
exact physical scattering amplitudes. 

The three-body relative Jacobi coordinates x,; and are 
used to describe the configuration of the three-body system in 
the cm. frame. So it is not required for proton to be infinitely 
heavy. The index i runs over all three arrangements of the sys- 
tem and corresponds to the spectator particle i. The positron 
Hydrogen Hamiltonian reads 

„ rj. , 6263 6361 6162 , . . ^ .-ir,. 
H^T-\ 1 1 \-igW2[X2). (10) 

Xi X2 X3 

Here T is the cm. kinetic energy, the index 1 corresponds 
to the positron, whereas 2 and 3 do to proton and electron, 
respectively, thus ci — e2 — e, 63 ~ — e. The last term igW2 
is the absorbing potential (|3j- For long-range potentials the 
modification |6] is needed, which consists in splitting of the 
Coulomb potentials into long-range Vl and short-range Vf 
pieces such that 

y/(x„y,) + l^/(x„y,) = y,^(x,) = (11) 

Xi 

This splitting is made in the three-body configuration space 
by a smooth splitting function Q i^i , Yi ) constructed in such 
a way that Q{K,,y,) = 1 if xJxq < {1 + yJyoY and 
0(x.i,yi) = Oifxi/xo > (l+yj/yo)" for some xo > 0, yo > 
and Q < v < 1/2. With such a Q the short- and long-range 
parts of the Coulomb potentials are defined as 

l^/(x„ y,) = Q{^,,y,)Vf'{^,y, = Vf - 

The Hamiltonian ilO\ is transformed then into 

iJ = i?' + y/ + V2' + igW2] =T+Vl+ + . 

Faddeev components of the wave function ^f"*" are defined by 
formulae 

= {E+ - HY^V{-^+ 

= {E+ - H^)-^{V2' +tgW2)'i>+, 

where = E + iO. It is straightforward to see that the sum 
of the components gives the wave function 
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and the components obey the set of MFE 



{E-H^ - K," - igW2)-i'^ = (Kj' + iffWa)^'^ (12) 

The important feature of equations il2\ . with regard to de- 
scription of the annihilation, is the fact that the two-body an- 
nihilation potential igW2 is incorporated into equations in two 
manners, i.e., in the diagonal part on the left hand side of the 
equation il2\ being responsible in the channel 2 for annihi- 
lation of Positronium after its formation, and in the coupling 
term on the right hand side supporting the annihilation process 
in the positron Hydrogen channel 1 through the coupling. 

The set of MFE (I12t treats the multichannel nature of 
the problem above the Positronium formation threshold cor- 
rectly. Namely, in the asymptotic part of the channel 1 where 
yi I oo and xi is confined within the Hydrogen atom vol- 
ume, the only nonvanishing potentials are Vj and Vf in the 
Vf = Vf and the first MFE takes the form 



combination 



0. 



{E -T -Vf)'^i 
So, it is obvious that asymptotically 

*+(xi,yi)^</>i(xi)e'^P-''\ 

where (j>i (xi ) is the Hydrogen wave function with the energy 
ei,sothati? = p^/2i/i + ei, and i/i is the e+-H reduced mass. 
In the asymptotic part of the channel 2, where |y2| oo 
and X2 is confined in the volume of the Positronium atom the 
nonvanishing potential is Vf + igW2 and the second MFE is 
reduced to 

{E-T-V^^- igW2)^t = 0- 

Hence, the asymptotics of the component can be written 

as 



*J(x2,y2) -^<^2(x2)e' 



E = Tpl/2v2 + e2 



where V2 is the p-Ps reduced mass, (j)2{y^2) is the Positronium 
wave function in presence of the absorbing potential, i.e., the 
solution to the equation 



A. 



X2 

2m 



e 

X2 



+ i.gM^2(x2)]'/'2(x2) = e202(x2). (13) 



Let us notice, that the first order approximation with respect 
to g gives for 62 the formula 



£2 



-M^°2\W2m 



with £2 and ip2 defined as the solution to (I13> if g = 0, i.e., the 
pure Coulomb problem for Positronium. 

The above arguments can be formalized rigorously with the 
use of integral equations. To this end it is convenient to intro- 
duce matrix notations 



H(5) = 









+ Vi + igW2 






Vf 


, w = 





0' 




■ 1 


0" 




, 1 = 




.^2 





W2 










1 



V = 

Thus, the MFE set takes the form 

[EI - nig) - V]* = igW^, 



(14) 



where = (5*]^,^^)^. The equation can be rewritten 

as 

* = *" + i.g[i;+I-H(c/) - V]"^W*. (15) 

Here '4'° is the solution to the homogeneous equation 

[£;i-H(.g) - V]*° = (16) 

which can be transformed to the integral form 

*" = * + [£;+I-H(5)]^iV*°. (17) 

The driving term $ obeys the homogeneous equation 

[£;i-H(.g)]* = 

and its components 4> = {^1,^2)'^ are given in terms of 
channel functions 4>i{xi) by formulae 

$,(x„yO =0,(x,)e'P'-'^'<5.fc 

which define through the index k = 1 , 2 of the Kronecker- 
delta the initial state of the e^—H system. 

From the analysis made in |5], |6] it follows that matrix in- 
tegral equations (I15> and ( I17> are well defined and possess the 
property that the driving term defines the solution uniquely. 
So that, the asymptotics of the solutions of il5\ and (I17> is 
entirely due to the asymptotics of the resolvent operators in- 
corporated in the kernels of these equations. The latter asymp- 
totics can be evaluated with the help of techniques developed 
in [5], 1 6], 1 1 1] which leads to the following result for the so- 
lution to (I15> as Ui ^ CO 



^'+(x,,y,,pfc) - (/),(x,)[e*P'-y-(5ifc + F,kiPryr,Pk)]- 

The amplitudes Fik have the structure which is similar to that 
of one channel case ^ 

The first term F^l is entirely due to the contribution of in 
( list and if i = k then F^j. is the elastic scattering amplitude 
whereas F^^. is the rearrangement amplitude if z 7^ fc. The 
second term is due to the absorbing potential. The direct an- 
nihilation amplitude in the e+— H channel is given by the 
formula 



(18) 
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Components ^'^"(p'j^) are solutions to ( I17> where the en- 
ergy should be taken on the lower rim to provide incoming 
boundary conditions, so that i?+ should be replaced by = 
E — iO, and the inhomogeneous term is set as $,j(xi,yi) = 
(/>,;(xi)e*P''^'(5,;i. The amplitude F^/. has the structure similar 
to that of ( I18> . It plays the role of the correction term either 
for the Positronium formation amplitude (fc = 1) or for the 
elastic amplitude of Positronium-proton scattering (fc = 2). 

Once amplitudes of all processes are determined, the cross 
sections of interests are defined by integrals 

aifc = ^ / dn = ^ / dn \F°, - ^gFl^f 

Vk J Vk J 

(19) 



TABLE II: Same as in TAB. I for L = 1 with same units. 



k 0.8000 0.8500 0.86315 0.8632 0.86325 

(721 0.4793 0.5621 0.1022[+1] 0.6847[-l] 0.4077 

(712 0.3633 0.3032 0.5167 0.3459 [-1] 0.2059 

crfi 0.581 [-6] 0.403[-6] 0.574[-6] 0.236[-6] 0.477[-6] 

al2 0.64[-6] 0.12[-5] 0.73[-6] 0.64[-7] 0.12[-6] 



'Ik 



Hi 

Vk 



dn [2g QmF^,Fl,+g^\Fl, 



(20) 



Here we have split the total cross section in the channel 1 
into two parts similarly to (0, so that cr*^. = aik + cr^f^, 
where ct°j. is the annihilation cross section. Formulae 
(El were used for calculations of the direct annihilation cross 
section af^ and the Positronium formation cross section a2i 
for the e+— H incoming channel, and the direct annihilation 
cross section (t°2 and the Hydrogen formation cross section 
ai2 when p— Ps is the incoming channel. During the calcu- 
lation we systematically used the smallness of the coupling 
constant g of the absorbing potential to simplify amplitudes, 
namely, from ( I15t it follows that '3/ = '3/° + 0(5) and from 
( fTTl it is immediately seen that \1>° = \I>°° + 0{g), where 
'3/°*' is the solution to ilH when g = 0. Hence, the ampli- 
tudes and cross sections in the leading order can be calculated 
with the solution of MFE for the pure three-body Coulomb 
problem. 

In tables below we list results of our calculations for (T12, 
a2i and cr"]^, <Ti2 cross sections for total angular momenta 
L = 0, L = 1 and L = 2 for five values of the energy of 
the spectator in the Ore gap. All calculations were performed 
on the base of a quintic spline algorithm for the solution of 
differential MFE equations from |12|. The numerical error 
for cross sections was estimated as not exceeding 5% for cr"i 
and 10%forcr?2- 

As a concluding remark we want to emphasis that the MFE 
formalism presented here is exact above the Positronium for- 
mation threshold, at the same time it is equivalent to the stan- 
dard one below the Positronium formation threshold iTsIl and 
it is also exact for the Positronium ion Ps- |9. 
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TABLE I: Ps formation cross section (721, H formation cross section 
(712 and direct annihilation cross sections an and (7i2 for L = 0. 
Momenta k are in units of 1/ao, and cr's are in units of vrag. 



0.8000 



0.8500 



0.8612 



0.8615 



0.8618 



(721 0.4965[-2] 0.5711[-2] 0.3239[-l] 0.826[-3] 0.3259[-2] 

(712 0.1137[-1] 0.9300[-2] 0.4982[-l] 0.1287[-2] 0.4972[-2] 

(7fi 0.299[-6] 0.329[-6] 0.399[-6] 0.332[-6] 0.336[-6] 

(7^2 0.30[-7] 0.37 [-7] 0.15 [-6] 0.53 [-8] 0.17 [-7] 



[1] A. Rich, Rev. Mod. Phys. 53, 127 (1981). 
[2] I. A. Ivanov, J. Mitroy, J. Phys. B: At. Mol. Opt. Phys. 33, L831 
(2000) 

[3] A. Igarashi, M. Kimura and I. Shimamura, Phys. Rev. Lett. 89, 
123201 (2002) 

[4] G. R Gribakin and J. Ludlow, Phys. Rev. Lett. 88, 163202 
(2002) 

[5] L. D. Faddeev, Sov. Phys. lETP 12, 1014 (1961); L. D. Faddeev 
and S. P. Merkuriev, Quantum Scattering Theory for Several 
Particle Systems (Kluver, Dordrech, 1993) 

[6] S. R Merkuriev, Ann. Phys. (NY) 130, 395 (1980); 

[7] P A. Eraser, Adv. At. Mol. Phys. 4, 63 (1968) 

[8] M. Charlton, J. W. Humberstone, Positron Physics (Cambridge 
University Press, Cambridge, England, 2001). 

[9] E. W. Schmid, H. Ziegelman, The quantum mechanical three- 
body problem (Vieweg, Braunschweig, 1974). 
[10] A. Igarashi, M. Kimura, I. Shimamura, N. Toshima, Phys. Rev. 

A 68, 042716 (2003) 
[11] S. L. Yakovlev, Theor. Math. Phys. 107, 513 (1996) 
[12] C-.Y. Hu, J. Phys. B: At. Mol. Opt. Phys. 32, 3077 (1999); 
A. A. Kvitsinsky, C-.Y. Hu, I. Phys. B: At. Mol. Opt. Phys. 
29, 2059 (1996); 

[13] A. A. Kvitsinsky, A. Wu, C-.Y. Hu, I. Phys. B: At. Mol. Opt. 

Phys. 28, 275 (1995). 
[14] A. A. Kvitsinsky, Yu. A. Kuperin, S. P. Merkuriev, A. K. Mo- 

tovilov, S. L. Yakovlev, Sov. J. Part. Nucl. 17, 113 (1986) 



TABLE III: Same as in TAB 1 for L — 2 with same units. 



k 


0.8000 


0.8500 


0.86549 


0.865505 


0.86552 


(721 


0.8711 0.1167[+1] 


0.8086 0.1788[-^1] 0.1306[+1] 


(712 


0.3961 


0.3779 


0.2425 


0.5362 


0.3917 


0^11 


0.283[-6] 


0.271[-6] 0.391[-6] 


0.271[-6] 


0.205 [-6] 




0.99[-6] 


0.12[-5] 


0.11 [-6] 


0.19[-5] 


0.58[-6] 



